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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM 1424 


ON THE THEORY OF ANISOTROPIC SHALLOW SHELLS* 
By S. A. Ambartsumyan 


1. INITIAL ASSUMPTIONS 

We shall consider a thin-walled, sufficiently shallow and aniso- 
tropic shell whose material, at each point, has a plane of elastic 
symmetry parallel to the middle surface of the shell. 

For the coordinate surface of this shell, we take the middle sur- 
face in curvilinear orthogonal coordinates a and B, coinciding with 
the lines of curvature. Let k-j_ = k-j_(a,B), k 2 = k 2 (a,,B) he the principal 

curvatures of the coordinate surface, and A = A(ct, B), B = B(a, B), he the 
coefficients of the first quadratic form. 

In regard to the shell, we make the following simplifying assump- 
tions : 

(1) The hypothesis of Kirchhoff-Love (ref. 2) shows that the recti- 
linear elements of the shell normal to the middle surface maintain their 
initial length after deformation of the shell, and remain rectilinear 
and normal to this surface. The error of this hypothesis, shown in ref- 
erence 3, has a value of the order of (8k) compared with unity, where 8 
is the constant thickness of the shell. 

(2) The parameters A(a, B) and B(a,B) are regarded as constants 
in differentiation (ref. 4). 

(3) Certain terms of secondary significance are neglected (ref. 5). 


"K teorii anizotropnykh pologikh oholochek." Prik. Mat. i Mekh., 
vol. XII, 1948, pp. 75-80. 

■^A solution of the analogous problem for a plate has been given by 
S. G. Lekhnitskii (ref. l) . 
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2 . EQUATIONS OP EQUILIBRIUM AND RELATIONS BETWEEN 
DEFORMATIONS AND STRESSES 

The conditions of equilibrium of an element of a shell, for our 
initial assumptions, are expressed by the equations 


B + A |§ + ABX = 0 
da o(3 

^9 Ss 

A w + B ^ + ABY = 0 




- (k^ + k 2 T 2 ) + 0%) + (AM 2 )] + Z = 0 


B 


bE 




( 2 . 1 ) 


£ “ A W ~ = ° 


dH 


3g 7 


- A 3F +B 33f + AB % = ° 


S-^ 4* Sg + ^2^2 = ^ 


J 


The last relation. In virtue of the formulas expressing the forces and 
moments through deformation of the middle surface, is an identity. 


For the deformations and the parameters of the changes in curvature, 
we have 


1 du 


_ _1_ d 2 v 

1 " “ A 2 3a 2 


6 2 


1 3v 
B 3p 


+ k 2 w 


03 


1 Su 1 Sv 

B 3]3 A 3a 


( 2 . 2 ) 


x 2 


1 d 2 w x _ _2_ d 2 w 

b 2 ^p2 AB 3adp 


(2.3) 


where u = u(a, p), v = v(a,p) are the displacements in the middle sur- 
face along the coordinate lines, and v = w(a,p) is the normal displace- 
ment. 


In equation (2.3) we neglected the components u and v in com- 
parison with the component w. Hence, these relations do not differ 
from the corresponding expressions for plates. This interpretation of 
the change in curvature for the general case was originally given by 
V. S. Vlasov (ref. 5). 
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Of the three differential relations given by A. L. Goldenveizer 
(ref. 6), for the case under consideration, only the last one is re- 
quired: 


1 


l b 2 < 


. „ , 12 1 d“CD ,1 

k 2Xl + k x x 2 + A 2 Sa 2 - AB + B 2 ap 2 


1 * 2 *1 

— - o 

O O w 


(2'.4) 


The equations of the generalized Hooke’s law in the chosen triortho- 
gonal curvilinear system of coordinates are 


a ” ~ A ll e a + A 12 e p + A 13 e y" + A 16 e ag 




= A 12 e a + A 22 e P + ^3^ + A 26 e a|3 
a T ~ A 13 e a + A 23 e P + A 33 e r + A 36 e aJ3 


T pr " A 44 e 3r + A 45 e ar 
T ar “ A 45 e |Br + A 55 e ar 


T ap A 16®a + A 26 e p + A 36 e r + A 66®a,p 


(2.5) 


In the case under considerations for a Y = 0, we have 


a a “ B ll e a + B 12 e p + B 16 e ap 
a p = B l2 e a + B 22 e p + B 26 e ap 
T aP = B 16®a + B 26®p + B 66®a,p 


> 


( 2 . 6 ) 


where, following S. G. Lekhnitskii 1 s theories (ref. l), there is intro- 
duced the notation 


B ik “ ^1^33 " A i3 A k3^/ A 33 = 


These stresses produce the following internal generalized forces: 
tangential T, S, bending, and torsional moments G, H, which, on the 
two principal sections a = constant and (3 = constant, have the form 



4 

Ti = 5 (Bh*i + Bi2 6 2 + Bigco) 

t 2 = 8(B 2 2 6 2 + B 12®1 + b 26“^ 

r3 

G 1 " " 12 ^ B ll x l + b 12*2 + B 16 T ) 

g 2 ~ “ jg ( B 22 x 2 + B 12 x 1 + B 26 t ) 

Si = - S 2 = S = 5( BggCD + + B 2 gS 2 ) 

' g3 T 

H l = ~ H 2 = H:= "i2 ( B 66 t + B 16 x 1 + B 26 x 2) ^ 

Substituting the values Gq, G 2 in equation (2.1), and H from 
equation (2.7), we obtain 

3 3 

% = - C(B ik )w N 2 = - B ( B ik> v ( 2 - 8 ) 

where 
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L (2.7) 




a 3 


1 a 3 

°< B lk) ’ B U i^ + 3B 16 a 2 b ^ • '-12 —66' ^2 ^2 ™26 j5 ^3 


+ B 

^ -Oort rr 


D(B ^ " B ^ & +3B ^ 55 ifc + (B ^> ^ SST : B “ 7 &J 


>(a) 


3. FUNDAMENTAL DIFFERENTIAL EQUATIONS 

For the unknowns, take u(a,p), v(a,P), and w(a,py. If 6g, 

co, k^, kg, and t from equation (2.2) and equation (2.3) are substituted 
in equation (2.7), T^, Tg, S, G^, Gg, and H can be determined as func- 
tions of u, v, and w. 

Further, substituting the values of the internal, forces in the equa- 
tions of equilibrium, and considering equation (2.8), there is obtained 
a complete system of equations for the three principal unknown param- 
eters, namely, u, v, and v. This system, following the work of V. Z. 
Vlasov (ref. 5), is presented in table I. 
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The equations (see 3.1, table I, connecting the unknowns u, v, w) , 
and boundary conditions make it possible to investigate the problem of 
the equilibrium of thin-walled, shallow, and anisotropic shells by the 
method of displacements. The integration of this system is, however, 
connected with very great difficulties. Making use of the method proposed 
by V. Z. Vlasov (ref. 5), for isotropic and anisotropic shells, the prob- 
lem can be reduced to a system of two simultaneous equations. 

We shall assume that X = Y = 0, that is, if the case of a surface 
with a normal load is considered. 

Setting 


*1 


JL_ <&p 
B 2 dp 2 


*2 


_1_ dfcp 

A 2 da, 2 


1 d 2 <P 
AB dad 3 


- (3.2) 


the first two equations of equilibrium are identically satisfied. Further, 
taking account of equations (2.3) and (3.2), we obtain, from equation (2.4) 
and the third equations of equations (2.l) 


^2 


1 d 2 w 1 d 2 w\ 1, . 

P 5P 1 P Sp) 8S 


■(’ 


(3.3) 


where 


ll(B “ ) " sp{ (BllBs6 - e ' 6) F + ps 5^- + 

^ B 11 B 22 “ B 12) - 2 ( B 12 B 66 " B 16 B 26)j ^“2 " ^ ' 2 + 

1 ^4 i ^4 'I 

2 ( b 22 B 16 " B 12 B 26^ ^3 ^^3 + ( B 22 B 66 “ ^4 ( 3 - 4 ) 


S2 = 


[(B^Bgg - B^g) (B 22 b 66 " B 26^ “ ( B 12 B 66 " B 16 B 26) 2 ] ( 3 ' 5 ) 

B 66 
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For the operators we introduce the notation 


« . n 2 A . is 2 

1 B 2 dp 2 * A 2 ba! 
Equation. (3.3 ) , in this case, assume the form 


= V V 
r v r v r 


(3.6) 


^l(®ik^ “ ^r w = 0 " " X 2 ^(®ik) w + ^ = 0 (3.7) 


t 

c\ 

c 


where tp » <p(a,p) is the stress function, analogous in the plane problem 
to the functions of Airy, and w = w(a, p) is the displacement function. 

From equations (3.7), for kq = k 2 = 0, we obtain the well-known 
equations for the plane stress state of a plate Lt(B -t> )<p =0, and for 
the bending of an anisotropic plate, L(Bjjg ) w = 12Z/6^. 

Thus, making use' of the mixed method of V. Z. Vlasov "(ref. 5), we 
obtain a more compact representation of the differential equations of the 
theory of anisotropic shells. The system (3.7) may be reduced to an equiv- 
alent single equation of the eighth order. We set 

w = L^B^)# <p = 62V r # (3.8) 

From the second of equations (3.7) , we obtain - 

+ p 2 4 * - P 2 (3.9) 


We note that this equation is a generalization of the equation given 
by V. Z. Vlazov (ref. 7), for isotropic cylindrical shells, and can be 
obtained by another method from the system 3.1 (table X), analogous to 
the method by which B. G. Galerkin (ref. 9), obtained the equation of the 
isotropic cylindrical shell (ref. 8). 


The internal forces, by equations (2.2), (2.3), (2.7), (2.8), and 
(3.8), are as follows: 


T x = 5S2 -L £-=• V r * 
1 B 2 dp 2 r 


12 ■ 7 h ^ 


S = - 82 B? SSfi 7 r* 


(3.10) 
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01 “ S ^ & + Bl2 + Bls ® ^] Ll(Bii)4 

°2 - n [ Bbb p & + ? & + »*• k &] (5 ' 1X) 

h - - § [fee h S55F + El6 ? f? + 526 ? 8?] 

N x = - ^ CCB^IrLCB^)* N 2 = - D (B i k) L l( B ik) # (3,12) 
For the displacement of a point of the middle surface, we have 

w - TaCB^)* (3.13) 


Sa 2 Sp 


/ -j ]_ 

u = - [(B U B 66 - B 2 6 ) k x + (B^Bge - B 16 B 26 )kgJ ^3 ^3 

[z(BiiB 26 - B^B^)^ - (B 22 B 16 - B 12 B 26 )k^J 

|[(BhB 22 - b| 2 ) - (B^Bge - B L6 B 26 )j k x - (B 22 B 66 - B§ 6 )k 2 | 

1 

( B 22 B 16 " B 12 B 26) k i ^3 ^ p 3 ) 

v = - g^-([(B 22 B 66 - B 26 )k 2 + (B 12 B 66 - B 16 B 26 )k^] p 
66 \ l— * 


1 


AB 2 dac>P 2 
(3.14) 


[ 2 ( B 22 B 16 - B 12 B 26 )B 2 ' (E 11 B 26 ’ B 12 B 1S ) ^l] ^ ^ 

|j(BllE22 - b| 2 ) - (B^Bee - B le B 26 )j kg - (BjjBee - b i 6 )*1^ 

1 

(Bn + B 26 - B 12 B 16 ) ^2 ^3 ^3 I 


’da 

1 S 3 # 


F + 


b&ar 

(3.15) 
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4. LOCAL STABILITY AND VIBRATIONS 


Repeat ing the considerations of v 7 in t , 

equations of* the local stabiH+tr ° f v - Z- Viasov (ref. 5 ) yields the 
form presented In table II. 7 M isotro P ic shallow shell in the 

to , C T£“ nt8 * “ d * are proportional 

( 3 . 7 ), we then obtain 7 ^ **» •*-**« 


^ L l( B ik)<P - V r w = 0 

** + “ L(BlkK ' fS £ i£ + 2S ° a sfe ♦ 4 $ §f] - 


( 4 . 2 ) 


equations 'for\hrstJess i StiJ ed r e Lrthe\r 1 r i ° n ° f tV ° simultaileo u 

This system can be reduced to a single . 5 , dis P la cement function w. 
a single function. From action ( 1 ^)°^“ ° f ^ 6lghth ° rder ' ** 


^l( B ik) B ( B ik)^ + ~§- Vy# 
b d r 


8 3 ^ t? Ic? + ^ 3® 33? + T 2 ^ ^Jb(Blh)* - 0 


( 4 . 3 ) 


tropic e hallow Shells an ? " /liration at anlso- 

computatioa, the Inertia forcea, ^%e 1 Sn“ the 

acceleratlon^of gravity* +f the , 8 h «U-< and g the 

respectively, ODtam, from equations ( 7 . 8 ) or ( 7 . 9 ), 

4 L l( B it )0 - V r w . 0 v r o * § L( Bli )w + Jj . Q (4 4) 

= 0 (4 .„ 

Form bfv.^vSoffre 6 ^: ™ «-t solved In this 

the equations ffiven +ho -fs ^ 1 2 > there is obtained, from 

Of anisotropic plates by S. GrSkSitSii^rSf ^8) ^ bllity and vibration 
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TABLE I* 


ufop) 

v(a,p) 

w(a,p) 

(3.1) 

, 1 a 2 *, 1 a 2 . B 1 a 2 

311 ? 5 ? +ZBlS AB 555P +B 68 -> ^2 

216 A 2 da 2 + +Be6> AB 5a5jf +B 26 32 Sp 2 

I^U k l + hs^ + |( B 16 k l + ® 2 S k 2 ^ 

X 

5 

^ ^ i? + f B L 2 + B 6 6 )ffi 5a5p +B 26 ^ ~2 

-1 ^2 i \2 -. n 2 

322 B 2 5P 2 + 2B26 AB 553p +B 66 *2 *a2 

f IP2&Z + B 12 k l^ + ^ B 26 k 2 + ®16 k l)^: 

Y 

8 

i(Bi 1 ki + B L 2 k 2 )^ ; + |<B 16 i 1 +B 2 6 l!: 2 )^ 

5 ( 222 ^ + %2 !l: l)^ + x( B 26 k 2 * B 16 k l)g5 

B 2 

C 8 !! 11 ! 2 + + ^22^2*) + 22 L ( B ik) 

Z 

"5 


«%k) - »U £ * * »L6 535- * +®6S) 


2 + 4 B 26 


1 a 4 

AB 3 dadp 3 


+ ®22 ~T ~^7 
B 4 5p 4 


I 


I 


■ _ I 1. 


1 • l. 


I 


*3Vt 5EL TOW 
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Table II. 


u(a,p) 

v(a,p) 

v(a,p) 

(4-1) 

B U ^ + 2816 A 353? + 

| - >2 
B 16 J2 ^2 + ( B 12 + B 66 )jS 353? + 

5(Bxi^i + * 1 ^ 1)55 + k l4 SS + 

b 66 Jj ^ 

1526 ^t?-* 1 * 260 

g( B 16 k l + B 26 k 2)gg + k l^ g gg 

^6^ ~g + ( B i2 + he)fe 335ff + 

Bg2 f&2 + A dW + 

g gg( B 22 k 2 + B 12 k l) + k 2 T 2 g ^g + 

k l k 2 B ° 

B 66 ^ - 44 

i ^( B 26 k 2 + B 16 k l) + k 2 S °I ^ 

j( B ll fe l + B 12 k 2^ + J ^ + 

g(B22 lE 2 + B 12 k l)^ + *2^2 g gg + 

j£ L (Blk) + ( B ll k l+^l^ k l k 2 + B 22 k 2) “ 

E^ 15 ! + B 26 b 2) ^ + *1^ i ^ 

i(B 26 b 2 + B-Lgk^ ^ + kgS 0 | ^ 

f® 7 8? + 3S3T +) * ^ ip) 
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